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1. Introduction 

The concept of 'open' quantum systems is a ubiquitous one in that all real systems 
of interest are open systems, each surrounded by its environment, which affects its 
dynamics. Caldeira and Leggett [1] used the influence functional approach developed 
by Feynman and Vernon [2] to discuss quantum dissipation via the paradigm of quantum 
Brownian motion (QBM) of a simple harmonic oscillator in a harmonic oscillator 
environment. The influence of the environment on the reduced dynamics of the system 
was quantified by the influence functional. Dissipation of the system originates from the 
transfer of energy from the system of interest to the 'large' environment. The energy, 
once transferred, is not given back to the system within any time of physical relevance. 
In the original model of the QBM, the system and its environment were taken to be 
initially uncorrelated. The treatment was extended to the physically reasonable initial 
condition of a mixed state of the system and its environment by Hakim and Ambegaokar 
[3], Smith and Caldeira Grabert, Schramm and Ingold [5], and Banerjee and Ghosh 
[6] among others. Haake and Reibold |7j , and Hu, Paz and Zhang [8] obtained an exact 
master equation for the quantum Brownian particle for a general spectral density of the 
environment. 

The spectacular progress in manipulation of quantum states of matter and 
applications in quantum information processing have resulted in a renewed demand for 
understanding and control of the environmental impact in such open quantum systems. 
For such systems, there exists an important class of energy-preserving measurements 
in which dephasing occurs without damping of the system. This may be achieved 
with a particular type of coupling between the system and its environment, viz., when 
the Hamiltonian H$ of the system commutes with the Hamiltonian Hsr describing 
the system-reservoir interaction, i.e., Hsr is a constant of motion generated by H$ 
[21 HO], [TTj . This condition describes a particular type of quantum nondemolition (QND) 
measurement scheme. 

In general, a class of observables that may be measured repeatedly with arbitrary 
precision, with the influence of the measurement apparatus on the system being confined 
strictly to the conjugate observables, is called QND or back-action evasive observables 
[T2l [131 [HI [15j [16]. Such a measurement scheme was originally introduced in the 
context of the detection of gravitational waves (T71 [18] . The dynamics of decoherence 
in continuous QND measurements applied to a system of two-level atom interacting 
with a stationary quantized electromagnetic field through a dispersive coupling has 
been studied by Onofrio and Viola [19]. In addition to its relevance in ultrasensitive 
measurements, a QND scheme provides a way to prepare quantum mechanical states 
which may otherwise be difficult to create, such as Fock states with a specific number 
of particles. It has been shown that the accuracy of atomic interferometry can be 
improved by using QND measurements of the atomic populations at the inputs to the 
interferometer [20J. QND systems have also been proposed for engineering quantum 
dynamical evolution of a system with the help of a quantum meter [21 J . We have 
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recently studied such QND open system Hamiltonians for two different models of the 
environment describable as baths of either oscillators or spins, and found an interesting 
connection between the energy-preserving QND Hamiltonians and the phase space area- 
preserving canonical transformations [22J. 

In this paper we wish to study the dynamics of decoherence in a generic open 
quantum system where the coupling between the system and its environment is of an 
energy-preserving QND type. The bath is taken to be initially in a squeezed thermal 
state, from which the common thermal bath results may be easily extracted by setting 
the squeezing parameters to zero. When the quantum fluctuations of the heat bath 
are squeezed, it has been shown by Kennedy and Walls [23J that the macroscopic 
superposition of states of light is preserved in the presence of dissipation. These authors 
have shown that the squeezed bath is more efficient than the thermal bath for optical 
quadrature-phase quantum measurements, and may also be used to prepare the states 
with low quantum noise in one quadrature phase, at least in the high-frequency regime. 
The advantage of using a squeezed thermal bath over an ordinary phase- insensitive 
thermal bath is that the decay rate of quantum coherences can be suppressed in a 
squeezed bath leading to preservation of nonclassical effects [24J. Such a bath has also 
been shown to modify the evolution of the geometric phase of two-level atomic systems 
[25] . In our present problem, we wish to systematically probe the effect of phase- 
sensitivity of the bath, and quantify the pattern of progressive decay of coherence of 
the system, both at high temperatures as well as arbitrary low temperatures, when a 
quantum nondemolition coupling is adopted. We wish to hence compare and contrast 
the quantum statistical mechanical features (viz. the nature of the noise channels) of 
the QND type of evolution with that of the dissipative evolution [HI [26], [27], [28], [29] of the 
general Lindblad form for a two-level system or the specific QBM form for a harmonic 
oscillator. 

The plan of the paper is as follows. In section 2 we obtain the master equation 
for a generic system interacting with its environment by a QND type of coupling. For 
simplicity, we take the system and its environment to be initially separable. In section 
2.1, the master equation is obtained for the case of a bosonic bath of harmonic oscillators. 
For the sake of completeness, we briefly compare the results for the oscillator bath 
with that for a bath of two-level systems in section 2.2. In section 3, we analyze the 
dynamics of decoherence, first by looking at the term causing decoherence in the system 
master equation for the bosonic bath of harmonic oscillators obtained in section 2.1, 
and explicitly solve it for the high-temperature and the zero-temperature cases. We 
then set up a quantitative 'measure of coherence' related to the linear entropy S(i), the 
dynamics of which is also studied for the zero- as well as the high-temperature cases 
for different degrees of squeezing of the bath. In section 4, the quantum statistical 
mechanical properties underlying the QND and dissipative processes are studied on 
a general footing for a two- level atomic system (section 4.1) and a harmonic oscillator 
system (section 4.2). For the two- level system, the dissipative process is taken to be that 
generated by a standard Lindblad equation while for the harmonic oscillator system, the 
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model studied is that of QBM. In section 5 we present our conclusions. 



2. Master equation 

Here we present the master equation for a system interacting with its environment by a 
coupling of the energy-preserving QND type where the environment is a bosonic bath 
of harmonic oscillators initially in a squeezed thermal state, decoupled from the system. 
We also take up the case where the environment is composed of a bath of two-level 
systems, and compare the two cases. 



2.1. Bath of harmonic oscillators 
We consider the Hamiltonian 

H = H s + Hr + H S r 



= H s + Y: to Ah + #sE 9k(h + b{) + H 2 S J2^- (i) 

k k k nuJk 

Here Hs, Hr and Hsr stand for the Hamiltonians of the system, reservoir and system- 
reservoir interaction, respectively. Hs is a generic system Hamiltonian which can be 
specified depending on the physical situation. b\, bk denote the creation and annihilation 
operators for the reservoir oscillator of frequency u>k, gk stands for the coupling constant 
(assumed real) for the interaction of the oscillator field with the system. The last 
term on the right-hand side of Eq. (1) is a renormalization inducing 'counter term'. 
Since [Hs, H S r] = 0, the Hamiltonian (1) is of QND type. The system-plus-reservoir 
composite is closed and hence obeys a unitary evolution given by 

pit) = e- iHt/h p(0)e iHt/h , (2) 

where 

p(Q)=p s (Q)p R (0), (3) 

i.e., we assume separable initial conditions. In order to obtain the reduced dynamics 
of the system alone, we trace over the reservoir variables. The matrix elements of the 
reduced density matrix in the system eigenbasis are 

-i(*£-i&)/ft£C92*/MO 



PUt) = e-^-W e 



x Tr 



R 



- iHnt/ VR(0)e iHmt / R J pl m (0), (4) 

where E n 's are the eigenvalues of the system Hamiltonian. In Eq. Pr(0) is the initial 
density matrix of the reservoir which we take to be a squeezed thermal bath given by 

p R (0) = S(r,$)p th S\r,<S>), (5) 

where 

Pth = II I 1 - e ^ hUlk ] e"^ 6 * 6 * (6) 
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is the density matrix of the thermal bath at temperature T, with (3 = 1/(&bT), 
being the Boltzmann constant, and 

6 t2 



S(r k , $ fc ) = exp 



(7) 



is the squeezing operator with r k , & k being the squeezing parameters [30]. In Eq. (j4]), 



H n = J2 \huj k b\b k + E n g k (b k + b\) 



(8) 



Following the steps of the derviation as shown in |Appendix A[ the reduced density 
matrix (4) of the system is obtained as 

Pnm(t) = e- l{E "- E - )t/h e 

_ 1) cosh(r fc ) + (e~ iuJkt - 1) sinh(r fc )e 2i * fc 2 



x exp 



Pnm(0). 



(9) 



Differentiating Eq. (Q with respect to time we obtain the master equation giving 
the system evolution under the influence of the environment as 



Pnm (^) 
where 

and 



h 



{E n - E m ) + if]{t){El - E 2 J - (E n - E m f^{t) 



Pnm 



v(t) 



E 



g 2 

h Ji 



1 



2 ^ n\ 



e iuJkt - 1) cosh(r fc ) + (e~ tuJkt - 1) sinh(r fc )e 



(10) 



(12) 



For the case of zero squeezing, r = $ = 0, and j(t) given by Eq. (fl2l) reduces to 
the expression obtained earlier [9j [TUJ [11] for the case of a thermal bath. It can be seen 
that r](t) fill I) is independent of the bath initial conditions and hence remains the same 
as for the thermal bath. Comparing the master equation obtained for the case of a QND 
coupling to the bath (JTDJ) with the master equation obtained in the case of QBM as in 
Refs. [261 123 [28] , where the master equation was obtained for the QBM of the system 
of a harmonic oscillator in a squeezed thermal bath, we find that the term responsible 
for decoherence in the QND case is given by 7(t). It is interesting to note that in 
contrast to the QBM case, here there is no term governing dissipation. Also missing are 
the various other diffusion terms, viz. those responsible for promoting diffusion in p 2 
and those responsible for diffusion in xp + px, the so-called anomalous diffusion terms. 
Also note that in the exponent of the third exponential on the right-hand side of Eq. 
([9]), responsible for the decay of coherences, the coefficient of j(t) is dependent on the 
eigenvalues E n of the 'conserved pointer observable' operator which in this case is the 
system Hamiltonian itself. This reiterates the observation that the decay of coherence 
in a system interacting with its bath via a QND interaction depends on the conserved 
pointer observable and the bath coupling parameters [10] • 
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2.2. Bath of two -lev el systems 

We briefly take up the case of a bath of two-level systems to illustrate in a transparent 
manner its difference with a bath of harmonic oscillators. The Hamiltonian considered 
is 

H = H s + H R + Hsr 

= H s + J2 UJ k&zk + H s J2 C kO- xk - (13) 

k k 

Since [H$, Hsr] = 0, the Hamiltonian ffT3l) is of a QND type. Starting from the unitary 
evolution of the entire closed system and then tracing over the bath variables, we obtain 
the reduced density matrix in the system eigenbasis as 

PnJt) = e-^-^^TrR [e^/V^VW] PLJO), (14) 

where 

H n = l^kVzk + E n C k a xk ] = O k {E n ). (15) 

k k 

Using the properties of the a z , a x matrices it can be seen that 

e iO k ( Em )t = CQS MEm )t) + + E m C k a xk ), (16) 

where 

u' k (E m ) = sjul + Elfll (17) 

Thus 

e iO k {E m )t e -iO k (E n )t = cos ^ Em )t) cos (uj' k (E n )t) 



sm(u' k {E m )t) sin {u' k {E n )t) , 2 2 
uj k (h m )uj k (E- ] 



sm(u}' k (E m )t)sm(u' k (E n )t) . . 

H . / Fi \ — / / i-i \ ^k^kK-^n — ^mj&zk&xk 

u k {E m )u' k {E n ) 
i cos (u' k (E m )t) sin KQgnjt) 

^k(En) 

icos(u] k (E n )t)sm(u k (E m )t) 
H 777; \ {u k or z k + &m<^kO~xk)- Uo) 

^fc(^m) 

Using (fT8|) in (jHj), it can be seen that only the first two terms on the right-hand side 
of Eq. (118j) contribute and the reduced density matrix of the system becomes 



Pnm{t) = e- i{En - Em)t/n Y[ 

k 



cos (uj' k (E m )t) cos (u' k (E n )t) 

PnrnW, (19) 



sin ( ^(E m )t) sin 

as also obtained by Shao et a/. [9J. We can see from Eq. (fl9|) that the reduced density 
matrix of the system is independent of the temperature and squeezing conditions of the 
reservoir, as may be expected from the structure of the Hamiltonian (fT3"|) . This brings 
out the intrinsic difference between a bosonic bath of harmonic oscillators and a bath 
of two-level systems. 
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3. Decoherence dynamics 

Here onwards we consider only the bosonic bath of harmonic oscillators in an initial 
squeezed thermal state. In this section, we analyze the decay of coherence of our generic 
system under various environmental conditions. We examine the decoherence term in 
the master equation for the reduced density matrix of the system. We also compute the 
dynamics of the 'measure of coherence' related to the linear entropy of the system for 
the zero- as well as the high-temperature cases for different degrees of squeezing of the 
bath. 



3.1. Dephasing in the system master equation 

In this subsection we examine in detail the term j(t) [Eq. (fl2l)]. This is the term whose 
time derivative is the decoherence-causing term as is evident from the master equation 
( TTUjl . To proceed, we assume a 'quasi-continuous' bath spectrum with spectral density 
J (a?) such that 

2 00 

Efl/M— > / dul{u)f(u), (20) 

and using an Ohmic spectral density 

I(u) = — ue- w/Uc , (21) 

7T 

where 7 and uj c are two bath parameters, we obtain 77 (t) in (ITT1) as 

rj(t) = tan" 1 ^). (22) 

7T 

In the limit u c t ^> 1, tan _1 (u; c t) — > | and 77 (t) — ► — y. Now we evaluate j(t) given 
in (fT2]) for the squeezed thermal bath for the cases of zero-T and high-T. 



T = : 

Using Eqs. (1201) . ( J2TT) in Eq. (TI21) and using the zero-T limit we obtain 7(7;) as 
, t s To _„x. /0 _m_/-, , ..2^ 7o „:_ U/0 _ N w [ (l + 4w c 2 (t-a) 2 ^ 



— cosh(2r) ln(l + cj 2 t 2 ) - — sinh(2r) x In 
2% An 



[1 + u 2 c (t - 2a) 



^sinh(2r)ln(l + 4aV 

4:71 



2 ^ 



(23) 



where t > 2a. Here we have taken, for simplicity, the squeezed bath parameters as 
cosh (2r(u;)) = cosh(2r), sinh (2r(a>)) = sinh(2r), 

$(w) = aw, (24) 

where a is a constant depending upon the squeezed bath. The decoherence-causing term 
7(7;) is obtained from the above equation as 



^-^cosh(2r) 

dt TT C ° Sh[2r) (1 + 



— sinh(2r) 

47T 



4.uj 2 c (t - 2a) 



1 + 4w c 2 (t - a) 2 1 + u 2 c (t - 2a) 5 



.(25) 
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Figure 1. ^ [ Eq . as a function of time t for different environmental conditions. Here 

70 = 0.1, u> c = 50, a = and temperature T = 0. The dashed and the solid curves correspond to the 
environmental squeezing parameter [Eq. i ]24[ )l r = and 0.4, respectively. 



We can see from the above equation that in the long time limit, 7(t) — > 
7 cosh(2r)/(7rt), and the terms proportional to the sine hyperbolic function, coming 
from the nonstationarity of the squeezed bath, are washed out. For the case of zero 
squeezing, we obtain from 



7(*) 



To 
2tt 



2.2\ 



ln(l + uj;t 



— x constant, 
2tt 



(26) 



because of the slow logarithmic behavior. 
As uj c — ► oo,7(t) in ff23l) tends to 

7(t) — ► — cosh(2r)A - — smh(2r)B, (27) 
2n Arc 

where 

A = lim ln(l + uj^t 2 ) = constant, because of the slow logarithmic behavior and 

l+4^( 4 _a) 2 )(l+4a 2 ^2)' 



UJ C >oo 



B = lim In 

U! c ^OO 

behavior. 



(l+^ (i _2a)2) 



constant, again because of the slow logarithmic 



High T : 

Using ([201) . fl2T|) in (fl"2l and using the high-T limit, we obtain 



<y k B T 
7(t) = — cosh(2r) 



irhuj t 
j k B T 



2a; c ttan 1 (d; c t)+ln 



1 



1 + tutt 



2+2 



sinh(2r) 



Au c {t — a) tan (2u c (t — a)) 



— Aoj c {t — 2a) tan 1 (cu c (t — 2a)) + 4aa; c tan 1 (2ouj c 



+ In 



'[1 +uo 2 c {t-2af 
[1 +Auj 2 (t-a) 



+ In 



Aa 2 uo 2 . 



(2? 
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where t > 2a. From (I28p we can obtain 7 (t) for high T and thermal bath with no 
squeezing, by setting r and a to zero, as 



lit) 



lok B T 

TfhuJr 



2u) c t tan 1 {oj c t) + In 



1 



l+uit 



2+2 



such that 



d 7 (t) 2 lo k B T 



tan 1 {uj c t). 



(29) 



(30) 



dt Tih 

This matches with the result obtained in Ref. [11] . For the case of a squeezed thermal 
bath, we can obtain j(t) from (1281) as 



dt 



cosh(2r) tan 1 (cu c t) ^° B — sinh(2r) 



— tan -1 (u c (t — 2a)) 



tan 1 (2u c (t — a)) 



(31) 



dy/dt 
0.2 
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Figure 2. ±g£L [Eq. J3TJ] as a function of time t for different environmental conditions. Here 
70 = 0.1, lu c = 50, a = and temperature T (in units where h = kg = 1) = 300. The dashed 
and the solid curves correspond to the environmental squeezing parameter [Eq. i l'24[ )l r = and 0.4, 
respectively. 

Figure 1 depicts the behavior of the decoherence-causing term, j(t) [Eq. §Z5§], for 
T = while Fig. 2 depicts its behavior for high-T [Eq. fl31~l) ]. with and without bath 
squeezing indicated by the parameter r. A comparison between the two clearly indicates 
the power-law behavior of the decay of coherences at T = and an exponential decay 
at high T. As u c — > oo, from ( 1281) we get 



7(0 



" ' ■ cosh(2r)t — 2 ^° B — sinh(2r)a. 



h 



h 



(32) 
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3.2. Evolution of the linear entropy - measure of coherence 



It is well-known that for a pure state p 2 {t) = p(t) and Tr[p 2 (t)] = 1. A mixed state 
instead is defined by the class of states which satisfies the inequality Tr[p 2 (t)] < 1. This 
leads to the definition of the linear entropy of a quantum state, S(t) = 1 — Tr[p 2 (t)], 
which is positive: S(t) > 0, and bounded: S(t) < 1. S(t) = for a pure state and 1 for 
a completely mixed state. We can set up a related 'measure of coherence' of the system 
following Ref. [9] as 



C(t) = Tr[p s (t)] 2 . 
If we assume the system to start from a pure state, 



m\ 



P s (0) 

then using ([9]) we have 

C(t) = Y,\Pn\ 2 \ Pm \ 2 e- 2 ^- E ^ t \ 

m,n 

where j(t) is as in ([12]) . 

The linear entropy S(t) can be computed from C(t) as: 



(33) 



(34) 



(35) 



S(t) = Tr p s (t) - (p s (t)) 2 



X 



C(t). (36) 

S(t) is plotted in Figs. 3 and 4 for a harmonic oscillator system starting out in a coher- 
ent state |a) [31], for temperatures T = and 300, respectively, and for various values 
of environmental squeezing parameter r. 



T = : 

Using Eq. fTl2|) in Eq. fl35l) and applying the T = limit, i.e., making use of (1231) . the 

measure of coherence is obtained as 

C(t) = £ \Pn\ 2 \Pm\\l + ^^ocosh^XM^V- 



70 sinh(2r)(£„-£ m ) 2 /(27r) 



;i+4^(t-q) 2 ) 
1 + u 2 c (t - 2a) 2 f 

70 sinh(2r)(_B„-_B m ) 2 /(27r) 



x ( 1 + 4a 2 o; 2 



For the case of zero-squeezing, r = = a, and C{t) given by Eq. (1371) becomes 



C(t) = J2\Pn\ 2 \p m \ 2 (uJ c t) 



-2 10 {E n -E m f /tt 



(37) 



(38) 



Here we have in addition imposed the condition uo c t ^> 1, which is a valid experi- 
mentally accessible domain of time. This agrees with the result obtained in [9]. It can 
be seen from Eqs. ( )37j) and (|38|) that coherences follow the 'power law' for T = 0. 



High T: 
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Figure 3. Linear entropy S(t) [Eq. d 36 D 1 as a function of time t for different environmental 
conditions: 70 = 0.1, U) = 1, w c = 50, a = 0, \a\ 2 = 5, and T = so that Eq. H23I I is used. The 
large-dashed, small-dashed and the solid curves correspond to the environmental squeezing parameter 
[Eq. r = 0, -0.3 and 0.4, respectively. 



Using Eq. f|T2|) in Eq. fl35|) and applying the high-T limit, i.e., using fl28|) . the measure 
of coherence is obtained as 

2 47ofc s T 



C(f) 



XI bn| 2 bm| 2 exp ^ -{E n - E n 

m,n ( 



7lk 



cosh(2r) tan 1 (u c t) — sinh(2r) tan 1 (2cu c (t — a)) 

Aa^ k B T 



+ sinh(2r) tan -1 (u c (t - 2a)) 



t — (E n — E„ 



7T 



Tl 



x sinh(2r) tan 1 (2oo c (t — a)) — 2 tan 1 (u c (t — 2a)) — tan 1 (2au; c ) 

x (1 + u ^y'yokBTcosh(2r)(E n -E rn ) 2 /(7Thuj c ) 

/M 2/ ,, _ . 2l 2\ 7ofcsTsinh(2r)(£: n -£; m ) 2 /(7r^ c ) 

/ [l + <(t-2a) 2 ] \ 
Ul+4a; c ^-a) 2 ]; 

X (1 + 4 :a 2 u j 2 ^-^okBTsmh(2r)(E n -E m ) 2 /{TThui c ) ^ ^39^ 

It can be seen from (l39l) that in the high-T case, the measure of coherence involves 
exponential as well as power-law terms. It is also evident that the terms dominating 
the temporal behavior of the coherence measure C(t) are 



bn| 2 bm| 2 exp I - (E n - E„ 

n,m ( 



cosh(2r) tan 1 {ui c t) 



sinh(2r)tan 1 (2u c (t - a)) + sinh(2r) tan 1 (u c (t - 2a)) 



t . 
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Figure 4. Linear entropy S(t) [Eq. d 36 D 1 as a function of time t for different environmental 
conditions: 70 = 0.1, £J = 1, ui c = 50, a = 0, |o| 2 = 5, and T (in units where h = kg = 1) = 
300 so that Eq. J28D is used. The solid, small-dashed and large-dashed curves correspond to the 
environmental squeezing parameter [Eq. i ]24[ )l r = 0, - 0.5 and 2, respectively. 

Thus in the high-T limit, the behavior of the coherences is predominantly exponential. 
In the long time limit (u c t —> oo), 

C(t) - £ |j9 n |> m | 2 exp \-{E n - E m ) 2 ^^- cosh(2r)t| . (40) 

n,m ^ J 

By comparing Figs. 3 and 4, it is evident that at T = (Fig. 3), the coherences 
stay for a longer time characterizing the power-law decay as opposed to the high-T case 
(Fig. 4), where the exponential decay causes the coherences to diminish over a much 
shorter period of time. Also evident is the effect of bath squeezing, characterized by the 
parameter r, on the coherences in the two temperature regimes. While in the zero-T 
case, the effect of squeezing remains over a longer period of time, in the high-T case 
it diminishes quickly. In this its behavior is similar to that of QBM of a harmonic 
oscillator system [271 EE] at high T. Another interesting feature that comes out is that 
in the zero-T regime (Fig. 3), by suitably adjusting the bath squeezing parameter r, 
the coherence in the system can be improved over the unsqueezed bath, as seen by 
comparing the small-dashed curve with the large-dashed one, representing the bath 
squeezing parameter ff2~4"l) r = -0.3 and 0, respectively This clearly brings out the utility 
of squeezing of the thermal bath. 

4. Comparison of the QND and non-QND evolutions and phase diffusion in 
QND 

In this section we make a comparison between the processes underlying the QND and 
non-QND (i.e., where [Hs,Hsr\ ^ 0) types of evolution for the system of a two- level 
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atom and a harmonic oscillator. We briefly consider the question of phase diffusion in 
the QND evolution of a harmonic oscillator. 

4-1. Two-level system 

Here we take the system to be a two-level atomic system, with the Hamiltonian 

Hs = (41) 

<y z being the usual Pauli matrix, to be substituted in Eq. (1). This is a common system, 
with a lot of recent applications, as for example, in the quantum computation models 

in ei Emu 



QND evolution 

In order to study the reduced density matrix of the system under a QND system- 
reservoir interaction, i.e., for using Eq. (|UJ), we need to identify an appropriate system 
eigenbasis. Here this is provided by the Wigner-Dicke states [351 EH EE] \j,m), which 
are the simultaneous eigenstates of the angular momentum operators J 2 and J z , and we 
have 

H s \j,m) = hum\j,m) 

= E j>m \j,m), (42) 

where —j<m<j. For the two- level system considered here, 3 = \ and hence 
m — — |, |. Using this in Eq. Q) and starting the system from the state 



|#))> = cos (!)|l> + e ^ sin (|W (43) 

the reduced density matrix of the system after time t is [25J 

s / cos 2 (f) |sin(^)e- i ^°)e-^) 2 ^)\ 

from which the Bloch vectors can be extracted to yield 

(a x {t)) = sin(^ o )cosM + 0o)e- (M27W , 
(o- y (t)) = sin(^ o )sin(cut + 0o)e- (M27(t) , 

(a z (t)) = cos(flo). (45) 

Here j(t) is as in Eqs. ff23|) . (1281) for zero and high T, respectively and o~ x , o~ y , a z are the 
standard Pauli matrices. It can be easily seen from the above Bloch vector equations 
that the QND evolution causes a coplanar, fixed by the polar angle #o> in-spiral towards 
the z-axis of the Bloch sphere. This is the characteristic of a phase-damping channel [38J. 

Non-QND evolution of the Lindblad form 

Next we study the reduced dynamics of the system (|4ip interacting with a squeezed 
thermal bath under a weak Born-Markov and rotating wave approximation. This im- 
plies that here the system interacts with its environment via a non-QND interaction 
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such that along with a loss in phase information, energy dissipation also takes place. 
The evolution has a Lindblad form which in the interaction picture is given by [311 [29] 

d 



dt 



p s (t) = 7o (iV + 1) (*- P s (t)* + - l -o- + a-p s (t) - ^(*K<7_) 

+ lo N (o-+p s (t)e_ - ^a_a + p s (t) - ^p s {t)a_a, 
- l0 Ma + p s {t)a + - 7o MV_p s (t)a_. 



Here 



N -- 
M 



and 



(46) 

(47) 
(48) 

e ftw/(fc fl T) _ 1 ' ( 49 ) 

where Nth is the Planck distribution giving the number of thermal photons at the 
frequency u, and r, $ are squeezing parameters of the bath. The case of a thermal bath 
without squeezing can be obtained from the above expressions by setting these squeezing 
parameters to zero. 70 is a constant typically denoting the system-environment coupling 
strength, and o~ + , a_ are the standard raising and lowering operators, respectively, given 

by 



N 



iV^cosh 2 r + sinh 2 r) + sinh 2 r, 



--sinh(2r)e i$ (2iV t/l + l) 



1 



tii 



|l)<0| = -(t7 JB +i < 7 y ) 



|o)(i| 



1 



o". 



(50) 



2 ri " y> ' 

In the above equation |1) is the upper state of the atom and |0) is the lower state. 
Evolving the system given by Hs from the initial state given in Eq. (j4"3"j) . using Eq. 
(I46| . we obtain the reduced density matrix of the system from which the Bloch vectors 
can be extracted to yield [25] 



Mt)) 



Mt)) 



,70 at 



1 (1 + cos$) 



-7o(2Ar+l+a)t/2 



<t*(0)> 



— sin $ sinh 



- 70 (2JV+l)t/2 



K(0)), 



1 + - (e 7oai - 1 ) ( I 



cos$) 



3 -7o(2AT+l+a)t/2 



K(0)> 



— sin $ sinh 

e - 7o (2iV + l)^ (0)) _ 



- 70 (2JV+l)t/2 



K(0)>, 



(2iv + r 



3 - 70 (2Af+l) 



(51) 



where 



sinh(2r)(2A^ + l). 



(52) 
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It can be seen from Eq. (1511) that the reduced density matrix p s (t) shrinks towards the 
asymptotic equilibrium state p aS ymp, given by 

Pasymp I g J i \pd) 

where p = \ 1 + ( 2 jv+i) ' ^ or ^ e case °^ zero sc L ueezm g an d zero temperature, this 
action corresponds to an amplitude-damping channel [38l [25] with the Bloch sphere 
shrinking to a point representing the state |0) (the south pole of the Bloch sphere) while 
for the case of finite T but zero squeezing, the above action corresponds to a generalized 
amplitude-damping channel [381 125] with the Bloch sphere shrinking to a point along 
the line joining the south pole to the center of the Bloch sphere. The center of the Bloch 
sphere is reached in the limit of infinite temperature. 

The above analysis brings out the point that while the case of the QND system- 
environment interaction corresponds to a phase-damping channel, the case where the 
evolution is non-QND, in particular where the evolution is generated by Eq. ( 1461) . having 
a Lindblad form, corresponds to a (generalized) amplitude-damping channel (for zero 
bath squeezing). This brings out in a very transparent manner the difference in the 
quantum statistical mechanics underlying the two processes. While in the case of QND 
interaction, the system tends (along the z-axis) towards a localized state, for the case of 
non-QND interaction, the system tends towards a unique asymptotic equilibrium state, 
which would be pure (for T = 0) or mixed (for T > 0). This can be seen from Figs. 5, 
where the effect of the environment on the initial Bloch sphere [Fig. 5(A)] is brought 
out. Figure 5(B) depicts the evolution under a QND system-environment interaction 
[Eqs. f l4"5l) ] while Figs. 5(C) and 5(D) depict the evolution under a dissipative system- 
environment interaction [Eqs. (1511) ]. While Fig. 5(B) clearly shows a tendency of 
localization along the z-axis, Figs. 5(C) and 5(D) illustrate the tendency of going 
towards a unique asymptotic fixed point. In Fig. 5(D), the presence of a finite $ ()48l) 
is manifested in the tilt in the figure. 



4-2. Harmonic oscillator system 

Next we take a system of harmonic oscillator with the Hamiltonian 

H s = hu (at a + ~) . (54) 

The Hamiltonian H s (1341 . substituted in Eq. (1), has been used by Turchette et al. [39] 
to describe an experimental study of the decoherence and decay of quantum states of 
a trapped atomic ion's harmonic motion interacting with an engineered 'phase reservoir'. 

QND evolution 

Noting that the master equation (ITOl) in the system space can be written equivalently 

as 

p s = - l -[H s , p s ] + ifi(t)[H 2 s , p s ] - 7(t) (H 2 s p s - 2H s p s H s + p s H 2 s ) , (55) 
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(B) 




Figure 5. Effect of QND and dissipative interactions on the Bloch sphere: (A) the full Bloch 
sphere; (B) the Bloch sphere after time t = 20, with 70 = 0.2, T = 0, £J = 1, UJ C = 40w and the 
environmental squeezing parameter [Eq. i ]24[ )l r = a = 0.5, evolved under a QND interaction [Eqs. 
145t ]: (C) and (D) the effect of the Born-Markov type of dissipative interaction [Eqs. H51ll | with 
70 = 0.6 and temperature T = 5, on the Bloch sphere — the x and y axes are interchanged to present 
the effect of squeezing more clearly. (C) corresponds to r = 0.4, $ = and t = 0.15 while (D) 
corresponds to r = 0.4, $ = 1.5 and t = 0.15. 



and substituting (|54|) in (|55l) we obtain the master equation for a harmonic oscillator 
coupled to a bosonic bath of harmonic oscillators by a QND type of coupling as 

p s = — itu^a, p s ] + ih 2 u 2 f](t) (a' a) 2 + a) a, p s 

- h 2 u 2 ^{t) \(a ] a) 2 p s - 2a ] ap s a)a + p s (a^a) 2 ] . (56) 

For clarity we transform the above equation into the form of a Q distribution function 
given by the prescription 
1 



Q{a) 



71 



-(a|p|a), 



(57) 
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where \a) is a coherent state. From the master equation (I5S1) . the equation for the Q 
distribution function becomes 

. d 



d_ 

ot 



Q 



d 



d 



— iuj a 



da* da 



Q + ih 2 LU 2 r](t) 



2(1 + aa*) 



a 



da" 



— a 



d_ 

da 



+ a 



a'' 



*2 



d 2 



da* 2 

d 



a 



d 



d 2 ' 

da 2 



da' 



+ a- — h a 

da 



*2 



Q - 
d 2 

da* 2 



h 2 umt) 



+ a' 



cP_ 

da< 



— 2aa* 



d 2 



dad a* 



Q. 



(58) 



Using polar coordinates a = £e , this equation can be transformed into 



d 2 



d9 " d£86 



Q 



d 2 



(59) 



Non-QND QBM 

To compare (|58|) with the equation obtained in the case of QBM of a system of harmonic 
oscillator interacting with a squeezed thermal bath, we start with the general QBM 
master equation [271 128] : 



d 

ih—p s (x,x',t) 



-h 1 



d 2 



d 2 



dx' 2 / 

„ (d_ _ d_ 

» dx dx' 
+ iDpp(t)(x — x') 2 p s (x, x' , t) 



2M \dx 2 

- ihr(t)(x - 



+ Y^mX 2 - X> 2 ) 



{^X j OC j ~fj j 



X 



p s (x, x', t) 



- h(D xp (t) + D px (t))(x - x' 




p s (x,x',t) 



- ih 2 D n 




p s (x,x' } t). 



(60) 



Here Y(t) is the term responsible for dissipation, D pp (t) for decoherence, D xx (t) promotes 
diffusion in p 2 , and D xp (t), D px (t) are responsible for promoting (anomalous) diffusion 
in xp + px. The details of these coefficients of the master equation (1601 can be found 
in [2"7l 12"%] . Here the coordinate representation of the density matrix has been used 
in contrast to the energy representation used in ffTUj) . Comparing fl60|) with ffTUj) we 
find that the QND coupling of the system with the environment makes the quantum 
statistical mechanics of the evolution much simpler. As already noted below Eq. ( |T2|) . a 
comparison between (jTOl) and (160]) shows that in the QND case there is a decoherence- 
governing term j(t), but no term responsible for dissipation. In contrast, the QBM 
case has dissipation and a number of diffusion channels as seen by the existence of the 
diffusion terms D xx (t), D xp {t) + D px (t) and D pp (t). 

Since Eq. (16"U|) is also obtained for a harmonic oscillator system (cf. f!51|) ). we 
proceed as before and obtain its corresponding Q equation as 



a 



_d_ 

da" 



— a 



d_ 

da 



Q 
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+ T r(t) 



(J 2 



da 



*2 



8F_ 
dot 



+ 2 



d 2 



dada* 



-2a 



d 



d 



da* 



-2a*— + 2a* 
da 



d 



da* 



+ 



ill 



muo 



D pp (t) 



d 2 



dada* 



1 

+ ■x 



d 

+ 2a— + 4 
da 

1 d 2 



d 2 



-^{D xp {t)+D px {t)) 



2 da* 2 
d 2 d 2 



*2 



ihmujD xx (t) 



d 2 



da 
1 

+ - 



da 2 



2 da 2 
Q 



Q 
Q 



d 2 



1 d 2 



Q. 



(61) 



dada* 2 da* 2 2 da 2 

From a comparison of Eq. (161j) with Eq. (15511 . it is evident that QBM is a more com- 
plicated process than QND evolution. Writing Eq. (1ST]) in polar coordinates does not 
simplify its structure, unlike the case of QND evolution where Eq. ( l59l) was obtained 
in a simple form in polar coordinates. This is a reflection of the fact that QBM is 
a more complicated process than QND as well as the fact that in the QND case, the 
master equation ( flOl) is obtained in the system energy basis which is more amenable 
to simplification in the Q representation (the Q-function being proportional to the di- 
agonal element of the density matrix in the coherent state basis) than the coordinate 
representation in which the QBM master equation (|60|) is obtained. 

Phase diffusion of the QND harmonic oscillator 

We now analyze Eq. (159|) to gain some insight into the process of phase diffusion in the 
case of a harmonic oscillator system coupled to its bath via a QND type of coupling. 
We take the long time limit. In this limit r)(t) — > (cf. remark below (122]) ). We solve 
for the Q distribution function for the zero and high temperature cases. 



T = : 

In the long time limit, rj(t) - 
to 

di Q = ^ 



0, and y(t) — > (cf. Fig. 1). Then Eq. §HD reduces 

(62) 



which has the solution 

Q(9,t) =e- xt e- X9/uJ , (63) 

where A is a constant. Equation (|62|) does not have the form of a standard diffusion 
equation in phase space - there is a drift term but no diffusion term. 

High T : 

In the long time limit, rj(t) — > and j(t) — ► joksT cosh(2r)//i (as can be inferred 
from Eq. (13T1) and Fig. 2). Then Eq. (|59l) becomes 

d _ d _ d 2 _ 
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with 

A\ = hu 2r y Q k B T cosh(2r). 
This has the elementary solution 



Q(0,t) =e~ at e~ M Cl e By + c 2 e- Bti , (65) 



,B0 i „ -Be 



where a, C\ y c 2 are constants, 
A : 



2Ai 
and 



uj AaAi 
2AV 

Equation fl6ll) has the form of a time-dependent diffusion on a circle. It does not have 
the form of a pure diffusion because of the presence of an additional Kerr-like term in 
the master equation ([SSD- A form similar to this arises in the phase diffusion model for 
the phase fluctuations of the laser field when the laser is operated far above threshold 
so that the amplitude fluctuations can be ignored [21] • Then the phase fluctuations due 
to random spontaneous emissions can be modelled as a one-dimensional random walk 
along the angular direction. In this sense it can be said that the QND Hamiltonian 
describes diffusion of the quantum phase [10] of the light field. From Eq. (1M|) it is 
evident that the diffusion coefficient is dependent on temperature T and the reservoir 
squeezing parameter r. In the high temperature and long time limits the dynamical 
behavior is that of a quantum mechanical system influenced by an environment that is 
modelled by a classical stochastic process, a situation that was studied in [TT]. That 
this is not so for the zero-T case suggests that a zero-T open quantum system cannot 
be simulated, even in the long time limit, by a classical stochastic bath. A detailed 
analysis of the phase diffusion pattern in QND types of evolution for the two-level atom 
as well as the harmonic oscillator system has been given in Ref. 



5. Conclusions 



In this paper we have studied the dynamics of a generic system under the influence of 
its environment where the coupling of the system to its environment is of the energy- 
preserving QND type. The bath is initially in a squeezed thermal state, decoupled from 
the system. We have compared the QND results with the case where the coupling is of 
a non-QND dissipative type for a system of two-level atom and a harmonic oscillator. 

For a bosonic bath of harmonic oscillators with a QND coupling (section 2.1), we 
have found that in the master equation of the system, though there is a term governing 
decoherence, there is no dissipation term, i.e., such systems undergo decoherence without 
dissipation of energy. For the case where there is no squeezing in the bath, our results 
reduce to those obtained in Refs. [9l [10l [11] for the case of a thermal bath. The reduced 
density matrix of the system interacting with a bath of two- level systems (section 2.2) 
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via a QND type of coupling is found to be independent of temperature [H] and squeezing 
conditions of the bath. This brings out an intrinsic difference between a bosonic bath 
of harmonic oscillators and a bath of two-level systems. 

We have analyzed the effect of the phase-sensitivity of the bath on the dynamics 
of decoherence, first by looking at the term causing decoherence in the system master 
equation for the bosonic bath of harmonic oscillators (obtained in section 2.1). We have 
evaluated the decoherence- causing term for the cases of zero and high temperatures, and 
also obtained its long time limit for both the cases. A study of the linear entropy S(t), 
which is an indicator of the coherences in the reduced density matrix of the system, 
clearly reveals (section 3.2) that in the high-T case, the effect of the squeezing in the 
bath is quickly washed out and the system loses coherence over a very short time scale. 
In contrast, in the zero-T case, the coherences are preserved over a longer period of time 
and the squeezing in the bath can actually be used to improve the coherence properties 
of the system. 

We have made a comparison between the quantum statistical mechanical processes 
of QND and non-QND types of system-environment interaction. For a two-level 
atomic system (section 4.1), it is seen that whereas the action of the QND system- 
environment interaction tends to localize the system along the z-axis indicative of, in 
the parlance of quantum information theory, a phase-damping channel [38], the non- 
QND interaction (epitomized by the Lindblad equation (146p ) tends to take the system 
towards a unique asymptotic fixed point, which for the case of zero bath squeezing would 
be indicative of the (generalized) amplitude-damping channel [38J. For a harmonic 
oscillator system (section 4.2), we have converted the master equation to the equation 
for the Q representation. This brings about in a very general manner the differences in 
the quantum statistical mechanical processes involved in QND and in QBM. The QBM 
process is much more involved than the QND one in that in addition to the decoherence 
and dissipation terms, it contains a number of other diffusion terms. In our analysis of 
the QND equation for the harmonic oscillator system, in the long time limit, we find 
a form similar to the one in the phase diffusion model for the fluctuations of the laser 
operated far above threshold when the amplitude fluctuations can be ignored. The phase 
fluctuations due to random spontaneous emissions can be modeled as a random walk 
along the angular direction. In this sense the QND Hamiltonian describes diffusion of the 
quantum phase of the light field. We find that while in the high-T case the situation can 
be modeled as a quantum mechanical system influenced by a classical stochastic process, 
it is not so for the zero-T case. The high-T Q equation resembles the equation of phase 
diffusion on a circle which would suggest a connection between quantum phase diffusion 
and QND evolution [H] . Our quantitative study provides a step towards understanding 
and control of the environmental impact in such open quantum systems. 
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Appendix A. Derivation of reduced density matrix (Eq. ((90) 



In this appendix, we present the steps leading to the derivation of the reduced density 
matrix p s nm (t) from (4) to ([9|). It follows from Eq. (jEJ) that 

e iH ( n k) t/h _ e -iE%gl(u k t-sm(u k t))/(h 2 ul) 



x D 



( E ^9k ( Wk t _ 1 

\ hut ^ 



iu k b]b k t 



where D{a) is the displacement operator, 



D(a) = e ab K- a * h 
and Y,H^ = H n . Similarly, 



(A.l) 
(A.2) 



= e iElal^ k t~sin^ k t))/(h^l) D f E n9k / 

V huoh ^ 



-iuJ k t 



x e 



-iu> k b\b k t 



(A.3) 

Now using Eqs. ( 1A.1I) and ( 1A.3I) in Eq. (BJ, and making use of the following properties 
of the displacement operator 

e ilJkb k bkt D(a) = D(ae iuJkt ) e iulkb * bk \ (A.4) 

D\a) = D(-a), (A.5) 

D\a)D(ae iu)kt ) = D (a(e iu)kt - 1)) e iaa ' sia{uJkt) , (A.6) 
the reduced density matrix in the system eigenbasis becomes 



PUt) = e-^- E ^ n e 



X 



nTr R [p R (0)D(fl k )]Pn m (0). 



Here 



9k = (E m — E n ) 



9k 



f^ k t _ ^ 



(A.7) 
(A.8) 



and Pr(0) is as in Eq. fl5j). 

The trace term in Eq. (1A.7I) is 

nTr R [p R (O)D(0 k )] = n Tr Rf S ( r k, < fk)pthS t (r k ,$ k )D(^ k ) 



k 

exp 



]J Tr R [p th D (6 k cosh(r k ) + 0* sinh(r k )e 2i * k 

/ phu k \ 

V 2 J 



1 k n ^k 



x 



[e iuJkt - 1) cosh(r fc ) + (e^ fct - 1) sinh(r fc )e 2i * fc 



(A.9) 
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Here we have used the following relation beween squeezing and displacement operators: 
S\r k , $ k )D(0 k )S{r k , $ fc ) = D (e k cosh(r fc ) + 6* k sinh(r fc )e 2i ^) . (A.10) 
Using Eq. I I A. 91) in Eq. ( 1A.7I) . the reduced density matrix p s nm {t) (EJ) is obtained. 
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